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We present numerical values for the self-energy shifts predicted by QED (Quantum Electrody- 
namics) for hydrogenlike ions (nuclear charge 60 < Z < 110) with an electron in an n — 3, 4 or 5 
level with high angular momentum (5/2 < j < 9/2). Applications include predictions of precision 
transition energies and studies of the outer-shell structure of atoms and ions. 



in 



in 
o 



> . 



I. INTRODUCTION 



The one-loop self-energy is the largest radiative correction in atoms and ions. It has been known for many years 
that for high nuclear charge Z, results obtained by a perturbation expansion in the number of interactions with 
the nucleus, i.e., in powers of Za are not accurate. There are many recent examples in which self-energy shifts for 
high principal quantum numbers and/or angular momenta are needed. For example, the outer-shell structure of very 
heavy elements is being studied to determine electron affinities or chemical properties. Similarly, transitions within 
the ground configuration of the Ti-like iso-electronic sequence (with outer-shell structure 3d ) have been measured 
for very high nuclear charges Z Jl|-0] , with results that differ systematically from theoretical calculations ||s|-|lT| that 
do not include the self-energy correction for the 3^/2 level beyond the Bethe logarithm and lowest-order electron 
■ anomalous moment. The same is true for systematic studies of magnesiumlike ions Yet it has been shown recently 
that effects beyond the lowest order are large at high Z, even for rather large n and angular momenta [l3] ]. Several 
calculations have thus been undertaken. Preliminary results have been reported for the 6s and 8^3/2 shells [ |l3) , 
' while Labzowski et al. [[uj and Yerokhin and Shabaev [jl5| have calculated a limited number of self-energy level shifts 
, with methods developed recently. In this paper, we provide accurate calculations for nuclear charge Z in the range 
00 ' 60 < Z < 110 for an electron in a level with n — 3, 4, or 5 with high angular momentum (5/2 < j < 9/2). The results 
increase considerably the number of atoms for which the exact one-loop self-energy is available. 

After the seminal works of Brown et al. [|l6| and of Desiderio and Johnson (jl7|, one of us (PJM) developed an 
efficient method for evaluating the self-energy level shifts in a framework where the electron-nucleus interaction is 
treated non-perturbatively. This method was first applied to the Is level of one-electron atoms [ p^|Jl9| . Later, the 
self-energy contribution to the Lamb shift of the 2s and 2^/2 states was studied j2(J, and subsequently, all levels 
with principal quantum number n = 2 were evaluated with higher accuracy |21| . Additional excited states with 
3 < n < 5 and angular momentum j < 3/2 were studied more recently in Ref. ]22| ], where the difficulties which arise 
for highly-excited states are described and solved for low angular momenta. Self-energy calculations for states with 
angular momentum j > 3/2 has been a long standing problem since the publication of the orig inal method ]l8| ]. 
In the present paper we extend the self-energy calculation methods developed in |18]-|22]] to arbitrary angular 
I- momenta. We derive general, numerically efficient formulas for the angular integrals that appear in self-energy 
expressions; these integrals were only known for states with j < 3/2. We perform calculations for 5/2 < j < 9/2 
and for principal quantum numbers n = 3, 4 and 5. The numerically efficient renormalization technique described in 
Refs. [ p3| , p4| is used, as its formulation is independent of the atomic level. This method yields precise results for the 
self-energy shifts; other groups have published numerical results based on various numerical strategies that give less 
precise values p5|~|27| . 

The implementation of our formulas on a computer is simple, since the only mathematical objects used in our 
final expressions are Bessel functions of the first kind and (squared) Wigner 3j symbols with all angular momentum 
projections set to zero. Our expressions are optimized for numerical calculations and have proved to yield very 
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accurate results. Our formulas for electrons with angular momentum j = 3/2 differ from known expressions 
because we have adapted our results to more accurate numerical calculations. 

The outline of the paper is as follows. In Sec. g we recall the principle of the self-energy calculation we use. In 



Sec. [II we summarize our formulas for the self-energy angular integrals, which constitute the main analytical result 
of the present work. In Sec. ^ we present numerical results for levels d 5 / 2 , f§/2, ./V/s; 27/2 an d 39/2 with principal 
quantum numbers 3 < n < 5. Detailed derivations of the formulas presented in Sec. Ill arc give n in Sec. v| while 



derivations by a different method, used for checking purposes, are presented in Sec. VI. Section VII concludes the 
paper. 

Throughout this article, we use the convention 53 i=Q = whenever a > b. 

II. SELF-ENERGY SHIFT FORMULA 

The expression for the self-energy shift of an electronic state n can be written for a large class of potentials V(x) 
as the sum 

£ SE = £ L + £ H (1) 
of a low-energy part £l and a high-energy part £n given (in units in which Ti — c — m e — 1) by 18 

£ L = -£ n - -P f £ " dz f dx 2 f dx 1 ^ n (x 2 )a l g(x 2: x l7 z)a m Vn (x 1 )(6 l ^V 2 ■ V x - V l 2 V?) ^F" ~ ( 2 ) 



& = dz J dx 2 J dx 1 (pl(x 2 )a^g(x 2 ,x 1 ,z)a^(p n (x 1 )^— 5m J dx ipl l (x)(3ip n (x), (3) 



and 



where b = —i [(£„ — z) 2 + ie] 1 2 , Re(6) > is the photon energy, and x 2 i = \ \x 2 ~ Xi\\ (e is an infinitesimal positive 
quantity). In these expressions, (p n and £ n are the eigenfunction and eigenvalue of the Dirac equation for the bound 
state n, and Q is the Dirac Green's function: Q{z) = (H — where 7i = a- p + V + P is the Dirac Hamiltonian. 

Indices I and m are summed from 1 to 3, and index \x is summed from to 3. The contour Ch extends from — ioo to 
— it and from + ie to +ioo. 

Since the high-energy part and the renormalization procedure described in Refs. pl]-p4| are already known for 
arbitrary angular momenta, this paper is concerned with the low-energy part (||). For a spherically symmetric 
potential V, a separation of the photon propagator and of Dirac wavefunctions into radial and angular parts yields 
the following expression for the low-energy part of the self-energy: 

(X (X r^ n /*oo poo 

£l = —£ n P / dz I dx 2 x\ I dx\x\ 

n ^ Jo Jo Jo 

2 

(x 2 , xi, z)fj(x 1 )A^ Kn (x 2 ,xi), (4) 

K i,j=l 

where the /i's are the radial components of the wavefunction ip n , 1 = 3 — i, are the radial components of the 
Green's functions, and K ^ are functions that contain angular integrations as well as the photon propagator. The 
Dirac angular quantum number of the electron for which the self-energy is calculated is denoted by n n . Detailed 
definitions of these notations can be found in Ref . |l8| . 

Expression (|]) can for instance be applied to a pure Coulomb potential V or to the self-energy screening correc- 
tion |2!| . The potential V affects the radial components / and G; on the contrary, the angular coefficients A are 
independent of V; this paper provides analytical formulas for these coefficients. 

III. FINAL FORMULAS FOR THE ANGULAR INTEGRALS 

In this section, we present our final formulas for the angular coefficients A Z J K introduced in Eq. (0) . Our results 



can be directly implemented on a computer. Derivations of the results presented here can be found in Sec. VA and 



Sec. VB 
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We present successively our results for the coefficients Aj. 1 K and A 12 K . The other two coefficients can be obtained 
through the following symmetries |f8| : 

A 21 = A 12 

A 22 = A 11 



A. Final result for Al] Kn 

1. Initial expression and notations 

We give in this section the numerically optimized form that we found for the angular coefficient Aj} K . The most 
general formula available for this coefficient is pl| 



Ai) Kn = \k\ j d£ 



-J-(1-?)P> h @P[ ®1±T0>) 

KK n P dp 

where P n is the Legendre polynomial of degree n, and where 

I k = \k+ 1/21-1/2 



(5) 



(6) 



denotes the orbital angular momentum associated with the Dirac angular quantum number k; the other functions are 
defined as: 



and 



p(xi,x 2 ,0 = \Jx\ + x\- 2zix 2 £ 



T(p) = sm(bp)/(b 2 p). 



(7) 



(8) 



The quantities n n , k, x\, x 2 and b are considered here as fixed parameters. Since distances x\ and x 2 often appear 
multiplied by the energy &, we define 

j/i = bx\ and y 2 = bx 2 . 
It is also convenient to introduce the notation 

J 2 EE 1(1 + 1) 

for the eigenvalues of the squared angular momentum operator in terms of the angular momentum I. 



2. Result of the analytical integration 
With the above notations, our result for the integration in Eq. (^J) reads: 



A% Kn = j2(u+m 



i=|i«-Z_„ 



2h i 1 q 1 Q n l Q ) ii{yi)ii(y2) 



KK n yi X 2 



(9) 



where the 2x3 matrices are Wigner 3j-coefficients, where the ji's are the spherical Bessel functions of the first 
kind pq], and where we have introduced a symbol for the following quantity: 
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C(a, fe, c) 



E 



i — \a — b\ 
i + | a — b\ even 

, otherwise. 







(2i + 1) ( J a 2 + J~ 5 2 - Jf) if a + 6 + c odd 



(10) 



Despite the non-symmetrical form ( fiol), the quantity £(a, 6, c) is totally symmetrical with respect to any permutation 



of its arguments, as we prove in Sec. V A 2 



3. Numerical implementation 



Equation (^) can be readily implemented on a computer. In Eq. (^), no term is singular (with respect to the 
parameters b, x\ and x 2 ). In fact, the expansion of the Bessel functions about the origin po| is such that ji(z) = 0(z l ); 
singularities could appear only through jo(z)/z, which is singular as z — > 0. But I = imposes that I + l K + l— Kn be 
even, so that jo(z)/z is actually never used in Eq. (^). 

Since the number of terms in Eq. (10) is smaller than c, we have decided to set c = l- K „ in Eq. (^J), because l- K „ is 
small, being the orbital momentum of the level for which the self-energy is calculated Jl8[ . Furthermore, we note that 
in Eq. (^|), for a given I, only one term inside the brackets contribute to A^} K because of parity selection on Wigner 
3j'-coefficients J|l| and in our quantity C(aJ), c) of Eq. (|l0|). 

We have also checked whether result (0) could yield strong numerical cancellations in the following ranges of 
parameters 



fee [0,1], Xl ,x 2 e [0,100], 
n I e [1,10] and |k| e [1,1000], 



(11) 



which are typical of the values used in numerical calculations |18[. We checked for possible cancellations with an 
arbitrary precision in the following two cases: (a) two numbers add up to a very small one (about 10~ 8 smaller than 
the two numbers); (b) a large number is added to a number which is much smaller (by a factor of 10 -9 ). Our checks 
have shown that no such numerical problems arise with the terms of Eq. (||) when the summation is done by starting 
from the higher limit l K + l- K „] this choice of summation order is numerically motivated by the over-exponential 
damping of the spherical Bessel functions of the first kind p3] : 



1 ( ez \ i 



(12) 



when I > ez/2; thus, summing the terms of the final result (|9|) with decreasing Vs allows one to include larger and 
larger terms in the total sum, which is necessary for obtaining accurate numerical sums. 

We have numerically tested formula (JsJ) against Mathematical against the alternate formula presented below 
[Eq. (|52|)] and against Mohr's implementation of his special cases k„| < 2; we have found excellent accuracy in 
the parameter ranges used in the low-energy part calculations of the self-energy [fl8| , ^i 22 1, i.e., in the range of 
parameters reached by our calculations of the self-energy in ions with Z > 10: 



fee [0,1], x u x 2 e [0,5000], 
|k„| e [1,5] and |k| e [1,5000]. 



(13) 



B. Final result for A^ Kn 
1. General integration result 



We give in this section a form of the angular coefficient A^ K which is optimized for numerical calculations. We 



start from the most general formula available |18| for this angular coefficient: 



A 



12 



1 



1 d 

pdp 



T{ P ) 



+ [x 2 Pi K .^) -^PlAm^Pi K JO-^iPi. n -AO}~~T(p) 

p dp p dp 



(14) 
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With the help of the notations defined previously, the (almost) final result of our calculation of Eq. (JT3) reads 



£ 

i=\i K - i Kn | 

l + l K + Ik.-, even 



3i{V\)ji{V2) N 
yi x 2 KK r , 



(2Z + 1)< 







UK J? 



jI-jL 



I — fa-, IfZn ) 



+ |«|(2f + 1) 







X 2 



' h) -j',ui^Jr • Jl JL)/K+3'i{yi)j'i(y2) x 2& 



Xl 



(15) 



which is the generalization to any atomic state of the angular coefficient Ajf K of the low-energy part of the self- 
energy Equation ( |l5| ) has a few crucial computational advantages: (a) all the sums contain a finite number of 
terms; (b) there are no large terms due to small b, x\ or x 2 ; in fact, it is easy to see that the only divergent quantities 
can come from ji=o(z)/z, but that this quantity is actually always multiplied by 0. We thus gain numerical accuracy 
compared to the published formulas pil, 22| for the special ±2. 



2. Numerical optimization 



As for A].\ , we have checked whether result ( [Tq ) yields strong numerical cancellations in the parameter ranges 
defined by ( pj|) . Our checks have shown that no such numerical problems arise with the terms of Eq. (filf) [evaluated 
in the order indicated by the parentheses, with summations done from the higher limit because of Eq. (12)], except 
sometimes for the minimum /, i.e., for 

^ ^min — \1>k ^K n I • 

However, we have cured this disease by noticing patterns in the way numerical cancellation appear: we have 



experimentally found four different areas of the (k, K n ) plane, defined by which of the terms of Eq. (15) cancel each 
other; the "terms" are defined here by separating terms with zero, one and two derivatives of Bessel functions: we select 
the four terms of Eq. (||) associated with \ji xnin (yi)/yi][ji mixx (y2)/x 2 ], j'i min (yi)\ji miri (y2)/x2\, b'i mi „(yi)/^i]ii min (2/2), and 
3L lB (Vi)fl min (V2)- _ 

By using explicit polynomial formulas for the orbital angular momenta l K and l Kn in each region, it is possible to 
obtain simplified expressions for the four terms of Eq. (|l5|); terms can then generally be grouped together [in a way 
that depends on the area in which (k, K n ) lies] through the Bessel identity [ p0| 

-3i(z) ~ fi(z) = 3i+i(*)- (16) 

We have found numerically that this identity can yield very strong numerical cancellations between the I = Z m i n terms 
ofEq. @. 

Thus, we have mathematically found the equations of the four cancellation areas, and we have used identity jl^ ) 
in order to express the final result in terms of the right-hand side of Eq. ( Jig ) instead of numerically calculating the 
left-hand side of the identity. The precise shapes of the areas in the (k, K„)plane are 



area 1 
area 2 
area 3 
area 4 



kk„ < and k(I k — l Kn ) < 
KK n < and k{I k — l Kn ) > 



KK n > and k(I k — l Kn ) < 
KK n > and k(I k — l Kn ) > 0. 

A numerical implementation of Eq. (|l5|) should thus calculate the terms corresponding to / = i m i n with the following 
expressions, which depend on the area in which (k, k u ) lies: 
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^k,k„ — 'mill) 



= < 



26| K |(2Un + i)r K " q k ^ in J iu te +i(yi)jw n +i(y 2 ) 

4 2 S JI = U) from Eq. (15) 

/ 1 I I ■ \ 2 
-|K|(2Z min + l) I q" q Q m 1 ji min +i(z/i) 

x {[juM/x^J^ J 2 + + 2bj' lm Jy 2 )} 

/ / / / ■ N 

-1 



\k\(21 x 







l(V2 



in area 1 
in area 2 

in area 3 
in area 4. 



(18) 



x {\juM/xi]{ J 2 _ + J 2 - JZJ/k + 2bf u j yi )} 



Furthermore, no factor of the form jo (z) / z, which diverges as z — > appears anymore in our formulas for A X2 K ^ ■ such 
a factor can only be found in the I = l m - ln term of result (fl5|); but Eq. (|l8|) must be used for this term. l mul = \l K — l Kn \ 
is obviously never encountered in the first two areas of Eq. (|l8|). In the last two areas of Eq. (|l8|), j (z)/z is used but 
multiplied by a factor zero, as can be seen by using the fact that kk„ > 0. 

In summary, Eqs. ( filf ) and ( |l8|) give numerically optimized formulas for the angle coefficient A 12 Kn of Eq. (|T^). We 
have numerically tested these formulas against Mathematica, against the alternate formula presented below [Eq. (|S0[)]. 
and against Molar's implementation of his special cases < 2; we have found excellent accuracy in the physical 
parameter ranges (O). 



IV. NUMERICAL RESULTS FOR THE SELF-ENERGY 

The formulas (^) , ( [l5| ) and ( presented above allow us to numerically evaluate the QED self-energy contribution 
£l formally expressed in Eq. (p"! The other self-energy contribution [£h_, Eq. (^)] can be calculated for an arbitrary 
electronic state by means of previously published methods pi] . We present in this section numerical evaluations of 
the self-energy of atomic electrons with principal quantum number 3 < n < 5 and with angular number 3 < \k < 5, 
in hydrogenlike ions with nuclear charge 60 < Z < 110. We used a Fortran implementation of formulas (^), (|l5) and 
( p^| ) , as well as previously published numerical procedures p9| , pl| |2jj2^,^,Q . 

Numerical results are most conveniently expressed in terms of the usual scaled self-energy F: 
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£ SE = -^-F(Zu)m c c 2 , (19) 
7r n 15 

where £se is the self-energy shift (|]). Our evaluations of the scaled value F are presented in Tables 0-|v|. When 
represented as a function of Za, our values display smooth curves, as can be seen in Figs. |l|-|^. 

From Figs. |l|-^, we notice that the scaled self-energy does not greatly depend on the atomic level n; put in other 
words, the n-dependence of the self-energy is quite well captured by the scaling 1/n 3 of Eq. (|l9|), as has been observed 
for s and d 3 / 2 levels [^3). On the practical side, the slow variation of the scaled self-energy F with respect to n allows 
one to use the values we present here as estimates for higher levels n. 

We have also grouped self-energy values by atomic level n in Figs. [| and |5[ We notice that states with identical 
orbital angular momenta have self-energy curves that look parallel on the scale we used. This property can be 
understood for low nuclear charge Z ', since states other than s states have the property that F is dominated in this 
region by a function of the form 

A 40 + (Za) 2 log[(Za)- 2 }A 61 

for Za 1, and where A^q and Aqi are the usual coefficients of the semi- analytic expansion of the self-energy 
shift |||. (The (Za) 2 A 60 contribution is negligible for Za -> 0: (Za) 2 = o[(Za) 2 \og([Za]- 2 )].) It turns out that 
for d,f,g... states, Aqi depends only on the orbital angular momentum (This is due to the smooth behavior of 
the wavefunction of such states at the origin j3(|.); as a consequence, the self-energy curves for states with identical 
quantum numbers n and I (angular momentum) are parallel in the limit Za — » 0. For the higher Z values of the 
results we present here, the expansion of F in Za is not supposed to hold; however our numerical values show that 
the difference between the scaled self-energies of two levels with the same n and I varies only on the level of a few 
percents in the range 60 < Z < 110. Precisely, the difference between the self-energy F for two states with identical 
n and I is well approximated by the value of this difference at the limit Za — > 0, which is known to be |3?J 

lhn Q [F n , K=Hl+1) (Za) - F n , K=l (Za)] = j^—^ , (20) 



G 



which approximates (to the level of a few percents) the splitting between states with identical n and I present in our 
results; it would be interesting to find an explanation of this property of the scaled self-energy F at high Z. At low Z, 
the splitting (^0|) is due to the anomalous magnetic moment of the electron |37j , and maybe does this effect dominate 
the splitting we observe between states of identical n and I in our high-Z results. 

Our results are coherent with the numerical results published by Yerokhin and Shabaev [M for nuclear charges 
Z = 74, 83 and 92, as shown in Tables 0and |V|. However, we note that our evaluations of the self-energy often lie below 



the values of ]2q] , and that our results have uncertainties smaller by about two orders of magnitude. Furthermore, 
the value that lies the furthest from their error bars is that of the 5rf 5 / 2 level for Z = 74, which is located relatively 
close ( — 1-2 standard deviations) from the previously published result pa]. 

V. CALCULATION OF THE ANGULAR COEFFICIENTS: METHOD FOR OPTIMIZED NUMERICAL 

EVALUATION 

Our analytical evaluations of A)^ R and A^ Kn Jl^l rest on a few cornerstones. 

First, the dependence on £ of the integrands of Eqs. (^) and ( jl4| ) is made simple by expanding the T function with 
partial waves. In fact, in Eq. (Q), p is the distance between the two interaction points of the self-energy, £ being 
the cosine of the angle between these two points ; the T function of Eq. (|J) can therefore be expanded in partial 
waves pcfl : 

T[p(x 1 ,x 2 ,0] = 1^(21 + l)m)Mvi)Mv*), ( 21 ) 



where we still have 



b 

l>0 



yi = bx\ and y 2 = bx 2 - 



Thus, integrals of terms of Eqs. ^ and ( fbi]) which do not contain derivatives of T are simply integrations of polynomials 
in £. And since d/dp can be calculated with d/dt; [see Eq. ( |2^ ) below], all terms are as well integrals of polynomials. 

Second, integrations of polynomials can be expanded in a simple form with the help of expansions on an orthogonal 
basis. Much of our work as thus been concentrated on obtaining formulas with polynomials in £ that can be simply 
expanded onto the basis of Legendre polynomials. 

Third, the sum over partial waves in Eq. (^lj) contains an infinite number of terms. Numerical calculations of 
angular coefficients Aj} K and ^ 2 Kn are more accurate if only a finite number of these terms contribute to the 
final result. This can be achieved by obtaining expressions in which we can apply triangular identities on Legendre 
polynomials (Legendre polynomials are a special kind of spherical harmonics); such a goal has driven our search of 
numerically efficient expressions for A RK and A^ Kn . 

Accordingly to these principles, expansions and scalar products in Eqs. (pij), (|30|), (|34|), and ( f45| ) are especially 
important to the following calculations. 

A. Integration in A\] Kn 
1. First steps 

We describe in this section our calculation (|^) of the integral over £ in the original expression of the angular 
coefficient A]^ Kn of Eq. (||) . 

Following Mohr J3f| , we evaluate the derivative 1/px d/dp of Aj} K in Eq. (||) by using a differentiation with respect 
to the £ parameter of the p function of Eq. (j^): in fact, for any function /, the definition of p implies that 



1 d tf \ 
pap 



= —4/^1,^,0]- (22) 



Starting from Eq. 
With Eqs. @ and' 



a first step consists in using the partial-wave expansion @ of T function defined in Eq. <M). 
22|), we thus obtain a form of A\ 1 K in which products of three Legendre polynomials (or 



derivatives of Legendre polynomials) must be integrated: 
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A% Kn = \K\Y,(2l + i)ji(yi)ji(y2) f <% 

l>0 



1 -1 



(i-e)PL(OPU (OPKO 



(23) 



where we have used the identity of orbital momenta l Kn —i — l- Kn - 

We can integrate over £ in the first term of expression (p3h with the well-known result H] 



d£P a (t)P b (0Pc(0 = 2 



a b c 




(24) 



2. Integral of derivatives of Legendre polynomials 
Integration of the second term of A}. 1 ^ in Eq. ( p3| ) is more difficult. We have found that 

r d£ (1 - emOPKOPiiO = C{a, b, c), (25) 



where £.(a,b,c) is defined in Eq. (10). The idea behind the derivation of this result was to consider the product 
(1 — £, 2 )Pa(t;)Pb(0 as t w0 coupled angular momentum eigenstates, so as to decompose this quantity onto the basis of 
Legendre polynomials, which we do now. The derivatives of Legendre polynomials are related to Legendre polynomials 
of orders — 1 and 1 by the following identities pa : 



v^em) = pho = -j-p^(o, (26) 

and the Legendre polynomials of a given order and degree are directly related to spherical harmonics : 



Y ltm (0,4>) = (-l) TO y (2 ^ ( ^ m) 7 )! j3 r(cosg)e^ ; (27) 

where the spherical harmonics Yi m (Q,(f)) is defined with the convention of Edmonds [ pi] , and where P" l {^) is the 
Legendre polynomial of degree I and order m |30| . We thus see that the left-hand side of Eq. (|2^) can be rewritten 
with spherical harmonics: in fact, we have from Eqs. ( p6| ) and (|2^) that 

(1 - cos 2 0)P' a (cos 0)Pl (cos 9) 
= -Jl Pi (cos 0)P^ (cos 0) 

= 4 V (2a4l + l) Fn ' l(M)n '- l(M); (28) 

the right-hand side of the above expression is constant with respect to cj> because the product of spherical harmonics 
has an orbital momentum projection of 0. The two spherical harmonics of Eq. ( p8|) can be coupled j3l]| : 



i>0 



Thus, the quantity (1 — ^ 2 )P' a {CjPb{Cj °f Eq. ( P5|) can be expanded over Legendre polynomials: with Eqs. (Pq), (|29| ) 
and (pTj), we obtain that 

(i-a^m'(e) = -y^E( 2i+1 )(i -1 o)(o o) p ^- ( 3 °) 



At this point, it is useful to have a closer look at parity selection rules in Eq. (|30j). Since only terms with a + b + i 
even contribute to the sum over i in Eq. ( |30| ) (because of the second 3j-coefficient), we can use the following formula 
in Eq. ©: 

a b i\_(abi \ J 2 — J 2 ~ 3 2 



1 -1 0) v° oy 2^ ' 
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provided that a + b + i is even |Q. On the other hand, if a + b + i is odd, then Eq. (|3l] ) does not hold; however, the 
second 3j'-coefficient of Eq. ( |30|) is zero in this case, so that we can safely plug Eq. (|31[) into Eq. (|30|): 



2 Wa 1 ^ b ' \ 



= E 2 -^(J*+J b 2 -J?)( a o o) 2 m). (32) 

z— |a — 6| 

This simple expansion over the Legendre polynomials is of particular interest in the sequel and will be used many 
times. 



3. Final steps 



With the help of Eq. (|32|), the integration of Eq. (|2jj) is straightforward if we know the coefficients of P{.(£) over 
the Legendre polynomials; they are easily obtained through an integration by parts: 



2 if i + c odd and i < c 
otherwise. 



(33) 



Taking into account Eqs. ( J32| ) and (33), we thus arrive at our final formula: 

c-1 



E 



i — \a — b\ 
i + | a — b\ even 

, otherwise. 



a b i 




(2i + l){ji + J h 2 -Jf ) if a + b + c odd 



(34) 



which is contained in Eq. (|Io| ) and (f25|). 

The angular coefficient 4^ obtained in Eq. (^3) can then directly be evaluated with the help of results ( p4] ) and 
([m|), and we directly obtain our final expression (||) for this coefficient. 

The evaluation of A^ K that we presented in this section also shows that only terms of Eq. (^) with I between 
\l K — l- Kn | and l K + l— Kn are non-zero; this fact, which was not obvious in the original formula, is explicitly expressed 
in our final result (ft). 



B. Integration in Al? Kn 

The angular coefficient A 12 Kn is more difficult to evaluate than A^} K , in particular because second-order derivatives 
of T appear in Eq. ( |l4|) . We show in this section how we calculate the integration over £ in Eq. (jlj) and obtain the 
final expression Eq. ( |15| ) [which must be combined with Eq. ( |l8| ) in numerical applications]. 

A first step consists in unifying the various terms of A 12 K in Eq. (|l4|); we can replace the Pi K = Pi_ K (£) of 

the last term of Eq. (|l|) by P lnn (£) and P^ (£) with the help of 

Pi-.J0 = -PLJ0+t p i.J0, (35) 
which can be easily deduced from pl[ , and which implies that 

x 2 Pi Kn (0-xiPl-«J0 

= ( X2 - Xl op lKn (o - ^(i - e)pi n (o- 

The complicated factor x 2 — x\^ can be removed with 

d 

x 2 ~x 1 ^p- — p, (36) 
ax 2 







which comes directly from the definition of p in Eq. By using 8/8x2 = {dp/dx^jd/dp, we finally obtain a special 
expression for the last term of A^ Kn in Eq. ( |l4| ) 



pdp 

Ox 2 K n Kn p dp 



(37) 



where P] K appears only in the left-hand side. Equation ( p7|) allows us to mix the various terms of the original 

A^k °f Eq. (U4) in a uniformed expression from which P^ _ t (£) has disappeared: 



(12 



N J d£ 



_9_ 

9x 2 



pdp 



KK„ ™ p Ctp 



1 d 



(38) 



It is useful to do a similar operation in which Pi K _ 1 (£,) is replaced by Pi„(£) and P/ (£), so that the terms are more 
uniform: by using again Eq. (B7D, but with k instead of K n , we obtain 



1 12 



i'.- / de 



_9_ 

9x 2 



*2-{i-e)piM)-^T- Pl M)^-T 

k P dp oxi 



-—(i-e)PLA0PL(0~T 

KK n ™ p dp 



-Xi 



+—(i-e)PL(o 



1 



1 d 



pdp 



9 1 d 



9xi p dp 



(39) 



As seen before, the derivatives d/dp are fruitfully calculated with Eq. (|22j), that we apply everywhere possible in 



9 
9x 2 

1 -1 



9 



(i-e^uiOPLjo^T 

KK n X1X2 Ot; 



K n X2 



(40) 



We note, however, that a second-order derivative 9 2 T/9£ 2 appears in the first term of the last line in Eq. (pl)|). In 
order to obtain a simple expression in which only first-order derivatives are present, we can use the following sort of 
integration by parts with Legendre polynomials: 



de(i-£ 2 )iY(Off(O/'(0 



dU(i + i)Pi(0g(0f(0 



d£(i-£ 2 )J?(0s'(0/(0, (41) 

L 

can easily be proved by integrating by parts and by using the differential equation of the Legendre polynomials p0[ |: 

d 



[(i-e)m)]=-ni+i)m)- 



(42) 



We thus transform the first term on the last line of Eq. ( flO| ) with Eq. ( f4l| ) [with I = l Kn , f = <9T/9£ and g 
(1 - £ 2 )P/. {£)]: this term contains 
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fx (i - ?) p L - £ 2 K (O^t = (i - ?)[JLPi*n (OK (0 + (0^»(0]^r, 



(43) 



where we have used the differential equation p^ ) on the g function, as well as the simple algebraic identity + 1) = 
K{K+l),i.e.,Jl=Jl 

We next use the partial-wave expansion ( pi|) of the T function in order to simplify integrations over angles; since T 
depends on £ only through Legendre polynomials [Eq. (pi])], inserting Eq. ( |43"| ) into the angular coefficient (40) yields 
only integrals of one of the following forms: 

dZP a (0Pb(0Pc(0 (44a) 

-l 

dzii-emoniOPcio (44b) 

i 

i 

2\ 



dca-n^m^K). (44c) 

L 

where the third Legendre polynomial of each line comes from Eq. ( pl| ) . We have already all the tools that allow us to 
calculate them: Eqs. (p|), (§51) and (^), along with the orthogonality relation fBOfl : 



(45) 



After using expansion (21) and the values for integrals (p|), we arrive directly at the final result ( p| ) for the coefficient 
A]^ Kn . We insist on the fact that the I = \l K — Z„ n | term of Eq. ( p| ) must, however, be evaluated with Eq. (|l8| ) in 
numerical calculations. 

As for -A^ 1 , the evaluation of A}? K _ made in this section shows that the sum in Eq. ( p| ) contains a finite number 
of terms: only the terms of Eq. ( p[ ) with I between \l K — l Kn \ and l K + l Kn are non-zero, as expressed in our final 
result (ph. 



VI. ALTERNATIVE METHOD FOR THE CALCULATION OF THE ANGULAR COEFFICIENTS 

In order to check our final results (||), ( p| ) and (|l8|), we derived and implemented independently a second set of 
expressions, that we compared to the results of the method presented in Sec. |v| over a wide range of arguments. This 
allowed us to limit the number of comparisons with Mathematical in fact, direct integrations of Eqs. @ and ( ph lead 
to very lengthy calculations. We found that the method presented in this subsection is accurate almost everywhere, 
but that it is slower than the method presented in Sec. and We present here the main steps of an alternate 
calculation of the integrals of Eqs. (](]) and (PI). 



A. Angular integration in A KtKn 

We evaluate Aj} K , starting again from the most detailed basic expression, Eq. (||). The idea of this second 
evaluation is to express the integrand solely in terms of Legendre polynomials of the variable f ; doing so will allow us 
to use the integration result of Eq. j2^). 

With this goal in mind, we remove factors 1 — £ 2 and £ in Eq. (^|) with the following Legendre polynomial identity 
[adapted to our particular form (|^) of orbital angular momenta] : 

-(i-£ 2 )it(0 = iU(0-tfU0. (46) 

and with the help of the Legendre recursion relation 

m(0 = ^ [Q + w+i(t) + iPi-im • (47) 

The T function of Eq. ( ^l|) can then be expanded with Legendre polynomials with Eqs. ( pl| ) and (p2|), and we can 
integrate once with Eq. (|24|); we thus get 
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<«„ = b\ K \j2(u + m (v2)ji (wo 



x <^ 2 







1 



K n y2yi J-i [ 



i 



2L + 1 

Formula ( [48| ) contains only one type of integral, namely 



(48) 



(49) 



The integrand of this formula can be transformed into a linear combination of products of three Legendre polynomials 
[easily integrated with Eq. (|24|)], by expanding the derivatives of Legendre polynomials over the orthogonal basis of 
Legendre polynomials: in fact, Eqs. (|45|) and (33) yield 



(50) 



where 



6 p (n, m) 



if n and m have the same parity, 

1 otherwise. 



In summary, we obtain an evaluation of the quantity T> m ' of Eq. (J49J) with the help of Eqs. ( p0| ) and fl24|): 

/l l—l n—l 

l — l n—l „1 

= E E( 2i + + Wp(n,j)6 p (l,i) / d(P m (0Pi(0Pj(0 



i=0 j=0 
l-l n-1 

2 EE( 2j + 1 )(^ + 1 ) ,5 p^' i ) ,5 p 

i=0 j=0 



m i j 




(51) 



With definition (fl9]), the final result of this section for A]} K ^ can be deduced from Eq. ( f48| ) and reads: 

< K „ = &NE( 2 ' + 1 ^ (yi) 



(>0 

^ 2 ' 'T 



n n y2yi 



V 



l:(i-«„> 



1 (/ 1)p l:(i_„,0 1 . p l:(/- K „,0 



(52) 



where the quantity P m ' is given explicitly by Eq. @. 

A numerical implementation of Eq. ( |5l| ) can make use of the fact that parity properties of the Legendre polynomials 
plf immediately impose 

D£(»>Q = if I + m + n odd. 

Moreover, numerical implementations are facilitated by the fact that the summation over I in Eq. ( j52] ) contains a 
finite number of non-zero terms — even though this is not obvious from the above form — , as noted for definition (|^) 
in Sec. [v| in Eq. (H|), the summation over I can be restricted to the range \l K — l- Kn \ ...l K — l— Rn . 
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B. Off-diagonal element 



To evaluate A 12 Kn we start from the most general published formula [Eq. (H)]. As for i^ 2 Kn , the idea of the 
derivation to follow consists in transforming integrations over £ into integrals of products of Legendre polynomials 
only, so that we can use the integration result (|2J). The first step is again accomplished by removing 1 — £ 2 factors 
with the help of Eqs. ([l6]) and (f47|), and by calculating derivatives of T with Eq. 



X2X1 J_i { Ot; K n " Ok, 



1 



i- K p_ K -i(t:)Pi n (o-^T(p) 



X2X1 



X2P-AO - xiPM 



X 2 P lnn (0-X!Pl_ Kn (0 



a 2 



(53) 



where we have used the simple identity on orbital momenta i K _i = 
A replacement of T by its partial-wave expansion (^lj) directly yields: 



.4 
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£(2/ + {y 1 ) Jl ( V2 ) [ ddp^iOPi^ (OPRO - -PA0PL (0*Y(0 

OX2X1 J —1 I K n 



1=0 
1 



(2/- 



1 



(/_ K + i)p i _ re+1 (e)P/ (0P/(0 



k„(2Z_ k + 1) 



(0P"(0 - Pi-A0P- Kn (OPl'iO PiA0Pi« n ®mt) + -PiA0P- Kn (0P"(0 

X\ " X2 

There are only two new categories of terms to evaluate. We thus define 

1.1 



V 



1-.1 



and 



V 2d _ 

m.n 



d^P m (0Pn^)P"(0- 



(54) 



(55) 



(56) 



Both these quantities can be expressed by means of an expansion of the Legendre polynomials derivatives with 
Eq. (p0|). We thus evaluate the first expression as 



V 



1-.1 



1-1 



I 

2^(2i + l)£ p (7,i 



-1 
i-i 



i=0 



to i n 




(57) 



where we have integrated products of three Legendre polynomials with (|2J. This formula is coherent with the fact 
that symmetry properties of Legendre polynomials |3l| ] show that expression ( |55| ) yields if to + n + I is even; this 
property should be used in numerical implementations. 

The function with a second-order derivative [Eq. (pq)] can be evaluated by using an integration by parts: 



V 2 ' 1 = 

m.n 



[P m (0Pn(0Pl(0}-l 



dz[p n (t:)Pm(z)] , pi(o 



[Pm(0Pn(0Pi(0ti- / dep»(e)j*(0iY(0 



d^P m (0Pn(0P!(0 



(58) 
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Using P m (l) = 1, P m {-1) = (-1)™ P/(l) = J(i + l)/2 and P/(-l) = + l)/2], we get 

,2:i ^ + 1) 



[l + ( _ ir +/+«]_^l:(-,0 + ^(«,0) 



(59) 



Symmetry properties of the Legendre polynomials and definition ( p6] ) show that this quantity is zero whenever l+m+n 
is odd, a fact that can be used in numerical calculations. 

Our final expression for the angular coefficient A 12 K thus reads: 



A 



^-£(2/ + l)j< (y 2 ) 



1=0 



1 



k„(2Z_ k + 1) 



I 



(2/— 



1) 2?. 



—V 2 ' 1 

x 2 



(60) 



where the various T> quantities are given explicitly in Eqs. (pl|), ( 571 ) and (|5^). We have checked the above formula by 
comparing Fortran and Mathematica outputs for many values relevant to the present work. Furthermore, numerical 
implementations should make use of the fact that it is sufficient to do the summation over I in Eq. ( |60|) only with / in 
the range \l K — l Kn | . . . (l K + l Kn ), because all other terms are zero [since Eq. ( |l4"| ) contains only these terms, as noted 
in Sec. Ml. 



VII. CONCLUSION 



We have obtained in Eqs. (0), (E5|) and (p"8|) analytic formulas for angular coefficients that appear in an efficient 



numerical method of calculation of the electron self-energy in hydrogenlike atoms 
angular momentum 1/2 < j < 3/2 were previously obtainable with this method |T~ 



18]; only results for electrons with 
231 ■ Recently developed numerical 



renormalization techniques |£3|,^4||29| allowed us to give numerical self-energy shifts for high-j states (5/2 < j < 9/2) 
of n = 3, 4 and 5 levels in the range 60 < Z < 110; our results are in agreement with recently available results [ |25| , 
but we provide values for many other nuclear charges Z , and the precision of our calculations is generally greater by 
about two orders of magnitude. These numerical results could for instance serve to include QED effects in many-body 
atomic calculations in atoms and ions with electrons of high angular momentum Q . 
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TABLE I. Scaled self-energy F(Za) for d 5/2 states (n = -3). 





n — 3 


71 — 4 




fin 




U.UO'iUoZl ^4 J 




00 




U.UODOO f 1 yQ) 




7(1 


n n^^i i (^ 






74 


0549734( A\ 






80 


0.0571916(5) 


0.0624150(5) 


0.0650174(4 


83 


0.0583498(4) 


0.0638430(4) 


0.0665677(4 


90 


0.0611703(5) 


0.0673377(4) 


0.0703682(4 


92 


0.0620040(5) 


0.0683750(5) 


0.0714981(9 


100 


0.0654457(5) 


0.0726763(5) 


0.0761897(5 


110 


0.0699331(6) 


0.0783217(5) 


0.0823592(6 



TABLE II. Scaled self-energy F(Za) for d 5/2 states (k = 3). 



z 


n = 4 


n = 5 


60 


-0.0203415(3) 


-0.0195626(3 


66 


-0.0201262(4) 


-0.0193062(4 


70 


-0.0199699(4) 


-0.0191179(4 


74 


-0.0198020(3) 


-0.0189139(3 


80 


-0.0195284(4) 


-0.0185761(3 


83 


-0.0193804(4) 


-0.0183914(4 


90 


-0.0190034(4) 


-0.0179156(3 


92 


-0.0188870(4) 


-0.0177669(4 


100 


-0.0183786(4) 


-0.0171113(4 


110 


-0.0176375(5) 


-0.0161389(5 



TABLE III. Scaled self-energy F(Za) for f 7/2 states (k = -4). 



z 


n = 4 


n = 5 


60 


0.0221590(3) 


0.0230942(5) 


66 


0.0225316(4) 


0.0235345(4) 


70 


0.0227973(3) 


0.0238498(3) 


74 


0.0230768(4) 


0.0241830(3) 


80 


0.0235224(4) 


0.0247170(3) 


83 


0.0237572(4) 


0.0249997(4) 


90 


0.0243367(4) 


0.0257006(3) 


92 


0.0245106(4) 


0.0259119(3) 


100 


0.0252427(4) 


0.0268060(4) 


110 


0.0262422(5) 


0.0280366(4) 
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TABLE IV. Scaled self-energy F(Za) for g 9 / 2 (« = —5) and g 7 / 2 (k = 4) states. 





AC = —5 


k = 4 


fin 






u u 


fl 01 98741 C ^ 




( U 


u.uizyooo^o ) 




74 


n m ^rizLi r/iS 


n m 940841^ 


80 


0.0131782(3) 


-0.0123450(4 


83 


0.0132500(3) 


-0.0123115(3 


90 


0.0134266(3) 


-0.0122290(3 


92 


0.0134794(4) 


-0.0122043(3 


100 


0.0137014(3) 


-0.0120991(3 


110 


0.0140037(3) 


-0.0119529(3 



TABLE V. Scaled self-energy F(Za) compared between previously published results and this 
work, for levels with principal quantum numbers n = 3 and 4. There is only one (small) incoherence 
between the results (3d 5 / 2 for Z — 74). 



74 
83 
92 



Ref. H 
This work 
Ref. H 
This work 
Ref. H 
This work 



n — 3 




n = 4 




3^5/2 


4d 5 /2 


4/s/2 


4/7/2 


0.0550(0) 


0.0598(4) 


-0.0198(4) 


0.0231(4) 


0.0549734(4) 


0.0596935(4) 


-0.0198020(3) 


0.0230768(4) 


0.0583(0) 


0.0639(3) 


-0.0194(3) 


0.0238(3) 


0.0583498(4) 


0.0638430(4) 


-0.0193804(4) 


0.0237572(4) 


0.0620(0) 


0.0684(2) 


-0.0189(2) 


0.0245(2) 


0.0620040(5) 


0.0683750(5) 


-0.0188870(4) 


0.0245106(4) 



TABLE VI. Scaled self-energy F(Za) compared between previously published results and this 
work, for levels with principal quantum number n — 5. There are only two (small) incoherences: 
5d 5 / 2 for Z = 83 and 5d 5/2 for Z = 74. 











n = 5 






z 




5d 5/2 


5/5/2 


5/7/2 


5^7/2 


5ff9/2 


74 


Ref. @ 


0.0628(6) 


-0.0184(7) 


0.0247(7) 


-0.0121(9) 


0.0135(9) 




This work 


0.0620673(3) 


-0.0189139(3) 


0.0241830(3) 


-0.0124084(3) 


0.0130416(3) 


83 


Ref. g 


0.0671(5) 


-0.0180(5) 


0.0254(6) 


-0.0121(8) 


0.0136(8) 




This work 


0.0665677(4) 


-0.0183914(4) 


0.0249997(4) 


-0.0123115(3) 


0.0132500(3) 


92 


Ref. §§ 


0.0719(5) 


-0.0175(5) 


0.0262(5) 


-0.0121(5) 


0.0136(6) 




This work 


0.0714981(9) 


-0.0177669(4) 


0.0259119(3) 


-0.0122043(3) 


0.0134794(4) 
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FIG. 2. Scaled self-energy F(Za) for / 5 / 2 states of levels n — 4 and 5. 
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FIG. 4. Scaled self-energy F(Za) for atomic level n = 4 and angular quantum numbers k = —3, 3 and -4. 
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FIG. 5. Scaled self-energy F(Za) for atomic level n = 5 and angular quantum numbers k such that |k| > 3. 
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